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ARTICLE INFORMATION ABSTRACT 

Article history: Stress might impair clinical performance in real life and in simulation-based education (SBE). Subjective 
Received 16 January 2017 or objective measures can be used to assess stress during SBE. This monocentric study aimed to eval- 
Received in revised form 18 June 2017 uate the effects of simulation of life-threatening events on measurements of various stress parameters 
Accepted 5 July 2017 


(psychological, biological, and electrophysiological parameters) in multidisciplinary teams (MDTs) dur- 
ing SBE. The effect of gender and status of participants on stress response was also investigated. Twelve 


Keywords: . emergency MDTs of 4 individuals were recruited for an immersive simulation session. Stress was assessed 
Simulation-based education orog ? coded ee 

Emercency by: (1) self-reported stress; (2) Holter analysis, including heart rate and heart rate variability in the tem- 
Multidisciplinary teams poral and spectral domain (autonomic nervous system); (3) salivary cortisol (hypothalamic pituitary 
Stress adrenal axis). Forty-eight participants (54.2% men, <7years of experience) were included. Measures were 


Scales performed at baseline (TO), after simulation (T1), after debriefing (T2), and 30 min after end of debriefing 
(T3). There was an increase in stress level at T1 (p< 0.001) and a decrease at T2 (p< 0.001). However, the 
variations of stress parameters induced by simulation (TO-T1 difference and T1-T2 difference) estimated 
by the three approaches were not correlated, while, as expected, Holter parameters were well-correlated 
to each other. Immersive SBE produced a change of stress level in all MDT members with no evidence for 
status effect but with gender difference. None developed a PTSD. These results support the hypothesis of 
a complementarity of the stress paths (collective reaction with increased stress level during simulation 
and a decrease during debriefing) but with relative independence of these paths (lack of correlation to 
each other). This study also suggests that because of the lack of correlation, stress response should be 
assessed by a combination of psychological, biological and electrophysiological parameters. 

© 2017 Australian College of Critical Care Nurses Ltd. Published by Elsevier Ltd. All rights reserved. 


1. Introduction 


Stress is defined as a state of real or perceived threat to 
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tively known as the stress response. Stress is both a psychological 
and a physiological phenomenon? and the physiological stress 
response to the stimuli is accompanied by an anxiety state (psy- 
chological stress). Excessive stress can lead to post-traumatic stress 
disorder (PTSD).? In health care activities workload, leadership 
issues, professional conflicts, and emotional care demands com- 
bine to generate stress.* This stress impacts performance, described 
by an inverted U-shaped curve in animal and human studies,>° 
and contributes to the definition of the arousal theory.’ Excessive 
stress and poor management can lead to human errors and decrease 
recognition of these errors in real life and in simulation-based 
education (SBE)® independently of professional experience.? The 
effect of stress on a medical team during emergency management 
of patients requires consideration, since it can compromise the 
safety of patients.!° Moreover, excessive stress interferes with the 
retrieval of conceptual knowledge and impedes problem solving. |! 
Therefore, evaluation of stress during SBE is of interest to enhance 
the pedagogical impact of SBE. While debriefing appears to be a 
crucial time in SBE,!*:!3 it is important to assess stress response 
during the simulation as well. Acute stress activates the autonomic 
nervous system (ANS), as can be apprehended by the accelera- 
tion of heart rate (HR) and modification in its variability (HRV) 
on Holter.!*!5 Stress also activates the Hypothalamic Pituitary 
Adrenal (HPA) axis activity, inducing biologic modifications such as 
an increase in cortisol that can be assessed by salivary cortisol (SC) 
in simulation.!° Finally, some events can induce perceived-stress; 
and a subjective approach to assessment of stress is commonly used 
in simulation to measure its impact.!’ However, the relationship 
between self-reported and physiological stress is not unequivocal? 
and may even be controversial.!® Consequently, stress response can 
be assessed by objective and/or subjective measures.!° A combi- 
nation of some biological, electrophysiological, and psychological 
measures of stress has been reported.2° However, most stud- 
ies have described the use of only one objective parameter of 
stress!®2!,22 or subjective parameter of stress.2°-2° Post-traumatic 
stress disorder (PTSD) can develop following a stressful event.’ 
Although SBE can produce excessive stress, the occurrence of PTSD, 
relating to experiences in the simulated learning environment, has 
never been studied. PTSD usually occurs from one week to one 
month after a psychologically traumatic event.2’ 

We hypothesized that all multidisciplinary team (MDT) mem- 
bers would experience stress during immersive SBE and that it 
would decline after debriefing. Because of their different physi- 
ological mechanisms, we did not expect a statistical correlation 
between subjective and objective stress measurements. We spec- 
ulated that stress response would differ according to gender and 
status within a team with a higher stress level in the leader. We 
assumed there would be no PTSD among the participants after a 
simulation session. 

The objective of this study was to evaluate the effects of sim- 
ulation of life-threatening events on measurements of various 
stress parameters (electrophysiological, biological, and psychologi- 
cal parameters) in MDTs. Secondary objectives were: (1) to evaluate 
the correlation between the different stress parameters (2) to study 
stress response according to status and gender; and (3) to search 
for the occurrence of PTSD. 


2. Methods 
2.1. Study 


This study took place in the ABS-Lab — Simulation Laboratory 
— INSERM (French national health and medical research institute) 
#1402, Faculty of Medicine of Poitiers, France. The study protocol, 
information form, and consent form were approved by the Comité 
de Protection des Personnes III de la region Ouest (Western France 


Person Protection Committee III) and were registered under the 
number 13.05.16. This study is the first part of a randomized single- 
center study on the relationship between stress and performance, 
registered by ClinicalTrials.gov under the number NCTO2424890. 
All details of study protocol were reported previously.”® The cal- 
culation of the number of subjects was based on an expected 
correlation of 0.50 between stress and performance. Four random 
samples of 12 participants for each of the 4 statuses present in 
an EMS team [Emergency Physician (EP), Residents (PGY i.e. Post 
Graduate Year), Registered Nurse (RN), and Ambulance Driver (AD)] 
were obtained. Each of the 12 MDTs was randomly constituted. 


2.2. Population 


All participants had less than 7 years of experience working in 
an EMS. This same level of experience allowed comparison in stress 
parameters between MDT members since stress depends on the 
competence and capacity of MDT leadership and co-workers.”° EPs 
worked in one of the hospitals in the Poitou-Charentes region (1.8 
10° inhabitants). They had obtained the University Diploma of Pedi- 
atric Emergency Procedures (University of Poitiers, France) after 
issuance of the American Heart Association and European Resus- 
citation Council guidelines in 2010.°° During the university course 
they had received identical training in pediatric insertion of intra- 
osseous (IO) access with performance assessment on the validated 
scale for simulated IO insertion.2! PGYs specialized in Emergency 
Medicine, were trained in pediatric emergency procedures: clinical 
training in a Pediatric Emergency Unit and/or the university course. 
All health care providers (RNs, ADs) from the EMS of the University 
Hospital of Poitiers had obtained the European Pediatric Immedi- 
ate Life Support degree over the two years preceding the session. 
All these EMS teams manage pediatric or adult patients in clinical 
practice. 


2.3. Intervention 


A high-fidelity manikin (SimNewB*, Laerdal® ) was used for 
immersive simulation of an infant with hypovolemic shock requir- 
ing IO access (Supplemental file 1). All sessions were scheduled in a 
simulated resuscitation room (Supplemental file 2) on the same day 
of the week at 2:00pm because of the cortisol circadian cycle. Each 
simulation — lasting 25-30 min — was preceded by briefing (15 min), 
and followed by debriefing (30-45 min). In addition, a 45-60 min 
“snack break” lasting until 5:00pm was planned after each simu- 
lation session to allow physiological variables to return to normal 
conditions. Stressful conditions were related to different sources:° 
the scenario itself, a realistic environment, and the untimely intru- 
sion of simulated parents in the resuscitation room. Four periods 
of 15 min were dedicated to saliva sampling and to filling out the 
STAI questionnaire. The saliva sampling was done at the same time 
for all participants. The chosen period for electrophysiological anal- 
ysis was 20 min long, starting from 10 minutes prior to the saliva 
sampling time, to 10 min after the saliva sampling time. 


2.4. Measurement of stress parameters 


2.4.1. Psychological parameters of acute stress and PTSD 

The STAI — State-Trait Anxiety Inventory scale?” is widely used 
in simulation.!® STAI included 20 questions with a 4-point Likert 
scale. Scores ranged from 20 to 80, with an acute anxiety response 
to stress for a score over 48 and anormal score around 34/80." Self- 
assessment of stress can also be evaluated by the Stress-O-Meter 
(SOM) scale, with a score ranging from 0 to 10.!74 The occurrence 
of PTSD can be assessed using the Impact of Event Scale-Revised 
(IES-R) on the 7th day after the event®° and the Post-traumatic 
Check-List Scale (PCLS) after one month.*° IES-R ranged from 0 to 
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and PCLS > 44/85.°° Participants filled out questionnaires in writing 
(E-mail). 


2.4.2. Holter parameters 

Holter parameters (HR (bpm), pNN50 and the LF/HF ratio) were 
obtained with the Synscope (Sorin Group’ ) software during 24h 
of recording, starting on the day prior to the simulation until the 
break after simulation. Cardiovascular activity was assessed by HR 
and HRV through time domain with pNN50 (proportion of succes- 
sive NN that differ by more than 50 ms divided by total number of 
NN)?” and the fast Fourier transformation (FFT) to perform spec- 
tral analysis and calculate LF/HF ratio (low and high frequency) 
that represents an estimation of “sympathovagal balance”.?®79 The 
sympathetic and parasympathetic activity always interfere with 
HR. HRV which reflect this interaction provide a measure of the 
sympathetic and parasympathetic activity induced by the stress 
response.*° It has been shown that acute stress, of which the exis- 
tence is confirmed by an increase in SC, involves a decrease in 
pNNS5O!° and an increase in LF/HF.*! 


2.4.3. Salivary cortisol 

SC is a faithful reflection of biological stress, with the major 
advantage that it does not require invasive sampling techniques.*2 
SC (wg/dl) was measured by an enzyme-linked immunosorbent 
assay (ELISA)*? with an ELISA kit (IBL international’, Hamburg, 
Germany).** Different parameters such as diet, nicotine and exer- 
cise may affect salivary cortisol but their impact was limited. The 
participants did not engage in any exercise the day before the sim- 
ulation. Moreover, as instructed by IBL procedure the participants 
did not eat, drink, chew gum, smoke, or brush their teeth during the 
30 min before sampling saliva, or rinse their mouths 5 min before 
sampling.*+ 


2.5. Statistical analysis 


Quantitative variables were described as mean and standard 
deviation (M+ SD). Categorical variables were summarized by fre- 
quencies and percentages. Stress scores were normalized on 100 
except for SOM, which was given on 10. The Kolmogorov-Smirnov 
test was used to check normal distribution for the measures 
assessed. Stress measurements were compared between base- 
line (TO) and the other times (T1, T2 and T3) using ANOVA for 
repeated measures (RM-ANOVA). Post hoc multiple comparisons 
among the groups were made using the Scheffe test. Regarding 
Stress-O-Meter comparison between T1 and T2 was performed 
using a pairwise student t test. Changes in stress level between 2 
times were assessed by calculation of the relative variation ((final 
score — initial score)/initial score). Correlations between the stress 
parameters were evaluated with Spearman correlation coeffi- 
cient. Stress responses were compared between socio-professional 
groups using a Kruskal-Wallis test and between males and females 
with a Mann-Whitney test. All tests were 2-sided, with type I error 
rate fixed at 0.05. All tests were performed using Statistical Analysis 
Software (SAS, version 9.3; SAS Institute, Cary, NC). 


3. Results 
3.1. Population 


Forty-eight participants (54.2% males) from the 12 MDTs who 
participated in high-fidelity simulation of emergency treatment of 
an infant with hypovolemic shock were considered in the analysis. 
The paramedics (RNs and ADs) were older than the EPs and the 
PGYs (3249 for EPs, 29+ 8 for PGYs, 42+ 13 for RNs, and 47+16 


Time variations in pNN5O (RV (TO-T1)/TO) 


Time variations in heart rate (RV (TO-T1)/TO) 


Fig. 1. Correlations between time-variations in heart rate and time-variations in 
heart rate variability estimated by time domain analysis (pNN50). 

pNN5S5O: percentage of RR’ intervals >500 ms; RV: Relative variation ((final score — 
initial score)/initial score); TO: baseline; T1: post-simulation. 


for ADs), and all of them had less than 7 years of EMS experience. 
The paramedics (RNs and ADs) had generally practiced in other 
departments before joining the EMS. 


3.2. Measurement of stress parameters (Table 1) 


Stress was studied on 12 MDTs of French EMS insofar as it 
might impair clinical performance in real life and in SBE. During 
the simulation session stress level increased with a significant vari- 
ation of stress measures between ‘pre-simulation time’ = baseline 
(TO) and ‘immediately after simulation time’ (T1): STAI, HR, LF/HF 
ratio, SC increased significantly between TO and T1 whereas pNN50 
decreased. ‘Immediately after debriefing time’ (T2) there was a 
significant decrease in the psychological and biological measures 
of stress, which all nonetheless remained above baseline (TO). HR 
was significantly lower and pNN5SO significantly higher 30 min after 
debriefing time (T3) than immediately after. The SOM score signif- 
icantly increased between TO and T1. Results are given in Table 1. 


3.3. Correlations between stress parameters (Table 2) 


Globally, there was no correlation between the different meth- 
ods regarding changes in stress parameters induced by simulation 
(TO-T1 difference, T1-T2 difference, and TO-T2 difference), with 
some subjects showing a more pronounced response in the psy- 
chological assessment, and others on the electrophysiological 
or biological parameters. The correlation between SC and STAI 
changes from TO to T1 (rho=0.33, p=0.02), and between SC and 
LF/HF ratio changes from TO to T2 (rho 0.41, p=0.006) were the 
only significant ones. Regarding the electrophysiological parame- 
ters, as expected there was a strong correlation between HR and 
HRV (Figs. 1 and 2), and also between time domain analyses and 
spectral domain analyses of HRV. 


3.4, Status and gender effects 


Changes in STAI scores, Holter parameters, and salivary cortisol 
measurements of stress induced by simulation were independent 
of status. While under baseline condition (TO) there was no gender 
difference in stress level, the stress responses induced by simula- 
tion between male and female were not similar. During simulation 
females displayed greater psychological stress with more recov- 
ery response after debriefing than males. STAI at T1 was higher in 
females (62.5 + 14.0 vs. 50.7 + 10.7, p=0.004). HR was significantly 
higher in females than in males at T1 (90.8+12.2 vs. 81.3410.9, 
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Table 1 
Objective and subjective stress parameters performed at TO, T1, T2, and T3, followed by between-time comparisons. 
Measures(mean + SD) TO T1 T2 T3 ANOVA Comparison Comparison Comparison 
F(3, 44) TO vs. T1 T1 vs. T2 T2 vs. T3 
SC 0.27 (0.14) 0.36 (0.24) 0.25 (0.15) ND 9.8 p=0.02 p=0.007 NA 
p=0.001 
HR 76.0 (9.9) 85.8 (12.4) 85.7 (11.7) 74.5 (8.6) 91.7 p<0.001 p=0.99 p<0.001 
p<0.001 
pNN50 16.3 (15.9) 9.4 (8.4) 10.1 (9.5) 15.8 (12.9) 15.9 p<0.001 p=0.98 p=0.01 
p<0.001 
LF/HF 4.5 (2.8) 6.3 (3.4) 6.1 (3.1) 5.93.7) 19.6 p<0.001 p=0.94 p=0.99 
p<0.001 
STAI 48.3(12.2) 56.4 (16.6) 40.6 (7.9) ND 35.8 p<0.001 p<0.001 NA 
p<0.001 
SOM ND 5.2 (2.1) 7.2 (1.8) ND NA NA t(47)=—5.98 NA 
p<0.001 


ANS: autonomic nervous system; HF: high frequency (Hz); HPA: hypothalamic-pituitary—adrenal axis; HR: heart rate (bpm); LF: low frequency (Hz); NA: not applicable; 
PNNSO: proportion of successive NN that differ by more than 50 ms divided by total number of NN; SC: salivary cortisol (g/dl); SOM: Stress-O-Meter (0-10); STAI: State-Trait 
Anxiety Inventory scale (over 100); TO: baseline; T1: post-simulation; T2: post-debriefing; T3: 30 min after debriefing time; ND: not evaluated; Values are means (standard 
deviations); Between-time comparisons were done by ANOVA for repeated measures and scheffé test except for SOM compared by t-test for paired measurements. 


Table 2 

Correlations between stress parameters. 
Correlation Between TO and T1 Between T1 and T2 Between TO and T2 

Rho p Rho p Rho p 

RV HR/RV pNN50 —0,64 <0,0001 —0,62 <0,0001 —0,74 <0,0001 
RV HR/RV LF/HF 0,01 0,96 0,20 0,17 0,38 0,01 
RV HR/RV SC —0,08 0,60 -0,13 0,37 0,01 0,93 
RV HR/RV STAI 0,25 0,09 —0,08 0,59 —0,03 0,86 
RV HR/RV SOM —0,15 0,30 NA NA NA NA 
RV pNN50/RV LF/HF —0,38 0,009 —0,56 0,0002 —0,50 0,0008 
RV pNN50/RV SC 0,06 0,67 0,03 0,85 —0,07 0,62 
RV pNN5O/RV STAI —0,05 0,75 —0,10 0,51 —0,01 0,92 
RV pNN50/RV SOM —0,01 0,93 NA NA NA NA 
RV LF/HF/RV SC 0,21 0,16 0,01 0,98 0,41 0,006 
RV LF/HF/RV STAI 0,10 0,51 0,11 0,47 0,06 0,38 
RV LF/HF/RV SOM —0,03 0,85 NA NA NA NA 
RV SC/VR STAI 0,33 0,02 0,29 0,05 —0,12 0,31 
RV SC/RV SOM —0,01 0,94 NA NA NA NA 
RV STAI/RV SOM 0,13 0,37 NA NA NA NA 


HF: high frequency (Hz); HR: heart rate (bpm); LF: low frequency (Hz); NA: not applicable; pNN50: proportion of successive NN that differ by more than 50 ms divided by 
total number of NN; RV: relative variation ((final score — initial score)/initial score); Rho: Spearman correlation coefficient; SC: salivary cortisol (jg/dl); SOM: Stress-O-Meter 
(0-10); STAI: State-Trait Anxiety Inventory scale (over 100); TO: baseline; T1: post-simulation; T2: post-debriefing. 
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Fig. 2. Correlations between time-variations in heart rate and time-variations in 
heart rate variability estimated by time domain analysis (pNN50). 

pNN5O: percentage of RR’ intervals >500 ms; RV: Relative variation ((final score — 
initial score)/initial score); TO: baseline; T2: post-debriefing. 


p=0.005) and T2 (89.4+9.7 vs. 82.34 12.5, p=0.03). LF/HF ratio 
was lower in females than in males (4.70+2.31 vs. 7.25+3.18, 
p=0.007). After debriefing STAI and LF/HF ratio decreased with 


more rapid recovery in females, whereas SC decreased with more 
rapid recovery in males. Relative variation of STAI between T1 
(post-simulation) and T2 (post-debriefing) was —0.18+0.18 for 
males vs. —0.32+ 0.19 for females (p=0.02). Relative variation of 
LF/HF ratio between T1 and T2 was 0.34 (SD=0.74) for males vs. 
—0.09 (SD = 0.45) for females (p = 0.03). Finally, relative variation of 
SC between T1 and T2 was —0.31+0.17 for males vs. —-0.21+0.15 
for females (p=0.03). 


3.5. Post-traumatic stress disorder 


No participant developed PTSD. Scores were transformed over 
100. IES-R score was 11.7+1.9 with a maximum at 34 and PCLS 
score was 19.0+0.8 with a maximum at 29. 


4. Discussion 
4.1. Main results 


The aim of this study was to evaluate the effects of simula- 
tion of management of life-threatening events on measurements 
of various stress parameters, pertaining to the impairment of per- 
formance by stress. In a simulation session with a model of infant 
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life-threatening emergency, stress level increased after simula- 
tion and decreased after debriefing for the 48 MDT care providers. 
STAI and SC measures significantly decreased immediately after 
debriefing, whereas the Holter parameters decreased later in the 
temporal domain but not in the spectral domain. There was no 
correlation between changes of the various stress parameters over 
time except for the Holter parameters between each other. There 
was no status effect. Females displayed a higher stress response 
during the simulation with a more rapid recovery for psychologi- 
cal and electrophysiological parameters. However, males exhibited 
a more rapid recovery in cortisol level. Although the simulations 
were stressful, there was no PTSD. To our knowledge, no previous 
study had searched for occurrence of PTSD during SBE. 


4.2. Primary objective 


The present study evidenced a rise in all parameters of stress 
during simulation of MDT management of an infant in shock. Simi- 
lar results were found in residents for perceived stress (Likert scale 
as SOM), HR and HRV, but not for blood cortisol after resuscita- 
tion simulation.*° The strength of the present study is to have 
studied stress response with multiple components as suggested by 
Laurent*° in high-fidelity simulation: perceived stress by 2 compo- 
nents (STAI and SOM) and objective measures of stress by Holter 
parameters and SC. Arora et al. suggested a similar approach to 
study stress response in the operating room. The authors pro- 
posed the use of the Imperial Stress Assessment Tool (ISAT), which 
is a combination of STAI, SC, and HR.*’ SC has the advantage of 
being a non-invasive method compared to the plasma cortisol 
determination used in other studies.*° Holter parameters included 
complete analysis of electrophysiological parameters (HR, HRV in 
time and frequency domains).!° Few studies have investigated 
these two axes in a clinical scene*® or in SBE.!® Furthermore, they 
were lacking in a standardized measure of stress.!*:2° In the lit- 
erature, usually when several measures of stress are employed 
they reflect only one pathway*® or the association of an objec- 
tive measure with self-report.*9 In SBE studies all stress pathways 
need to be assessed because variation in the systems could be 
different.°° Moreover, it would be interesting to estimate the task 
level for participants and thereby predict excessive stress during 
training.°'! The present study showed an activation of all stress 
axes during the simulation session, with a decline occurring after 
debriefing. Decrease in stress level after debriefing has never been 
described. This could be related to the good-judgement debrief- 
ing procedure per se.°* But in the absence of a control group 
(without debriefing), spontaneous decrease in stress measures 
after simulation cannot conclusively be ruled out. After debrief- 
ing, the decrease in SC and perceived stress was faster than for 
Holter parameters. Moreover, while the spectral domain param- 
eters remained stable, there was a decrease in stress parameters 
calculated in the temporal domain. These results suggest a more 
complex relationship between stress, HR and its different estimated 
variability parameters (temporal and frequency/spectral), and the 
evidence of LF/HF variation such as sympathetic modulation®* and 
parasympathetic system withdrawal.°* Nevertheless, the present 
results reinforce the idea that scheme of frequency bands can- 
not fully represent a model of ANS response to stressors,?° as 
multiple factors are liable to interfere with sympathovagal bal- 
ance. 


4.3. Secondary objectives 


Despite all stress markers increasing during simulation and 
declining after debriefing, there was no correlation (except for HR 
and HRV) between their absolute values and their variations. How- 
ever, for ANS parameters there was a higher correlation between 


HR and pNN50 (i.e. temporal domain analysis) than HR and LF/HF 
ratio (i.e. spectral domain analysis). Our results suggested that the 
use of LF/HF ratio as a surrogate for sympathovagal balance was 
contentious. The LF/HF ratio may not accurately measure cardiac 
sympathovagal balance, particularly because of the poor relation- 
ship between LF power and sympathetic nerve activation described 
in the litterature.°°°” Surprisingly, Hunziker found a correlation 
between plasma cortisol, HR, HRV, and perceived stress on base- 
line that contrasted with the absence of correlation after simulation 
and debriefing except for HR and HRV.*° Variations of objective 
measures of stress differed between each other in two dimensions: 
cortisol rate of change differed from that of electrophysiological 
Measures and variation over time was not similar. Higher per- 
ceived stress during simulation was not correlated with ‘biologic 
stress parameters’. This study seemed to evidence that the link 
between stress pathways is not so simple. On the contrary, the 
mechanisms of biological — objective - and subjective stress param- 
eters expression are complex, as is suggested by contradictory 
results,?:!6!9.23.33 Absence of correlation might be explained by 
differences between stress response circuits. Stress response is 
regulated by neural circuits impinging on the Hypothalamic Pitu- 
itary Adrenal (HPA) stress axis! and the sympathetic branch of 
the ANS.°® ANS leads to a rapid response for coping mechanisms 
by modulation of electrophysiological parameters and released 
catecholamines whereas cortisol is issued from a slower cascade 
activated by HPA.*° Cortisol levels peak at 30-40 min after the 
onset of stressors.°° Indeed, even if HPA and ANS are interrelated 
and interact,°° there was no correlation in spite of a non-equivocal 
increase in all measures of stress during simulation. Some authors 
have described ANS parameters as reflecting an ‘effort or challenge’ 
(alertness and action proneness) component of stress and HPA 
cascade and a ‘distress’ (negative feeling) component.*° One expla- 
nation of the absence of observed correlation between objective 
and subjective measures of stress parameters could be the inter- 
play of the neurohormonal axis, circadian rhythms, and individual 
variability in chronotropic response to stress.°! Moreover, Metzen- 
thin et al. found that concurrently self-reported stress like SOM was 
related to SC but that this was not the case if it was retrospectively 
determined.*9 In short, face with a lack of correlation between 
stress parameters it may be advisable to investigate parameters 
off all pathways, and also to pay particular attention to subjec- 
tive parameters when investigating stress response in simulation 
courses. 

The present study found that as previously described,!*:°2 
all participants exhibited a similar stress response irrespective 
of their status. This suggests that stress might be assessed in 
all team members and rather than focusing on one of them 
(leader). Indeed, excessive stress in one team member may impair 
teamwork.°? As expected, there were gender differences according 
to stress response. However, females expressed a higher auto- 
nomic response on baseline than males, whereas usually, stress is 
higher in males.°*°° But other authors have suggested that there 
is no gender difference or that stress response could be greater 
in females.°°°° Similarly and as previously suggested,°” females 
expressed more anxiety during simulation. Analysis of relative vari- 
ation of SC indicated that SC decreased more quickly among males 
after debriefing. McGraw et al. suggested larger changes in SC level 
in males and more rapid recovery after simulation.°® However, 
to our knowledge, this variation after debriefing has never been 
described. Conversely, decrease of STAI after debriefing was more 
pronounced in females than in males. Finally, relative variation 
of LF/HF ratio varied in the opposite way in females and males 
reinforcing the hypothesis that multiple factors might interfere 
with sympathovagal balance and gender difference. All in all, the 
heterogeneity of these results highlighted the complex nature of 
gender-specific stress response.°? Our results showed that gender 
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could be treated as an independent variable in all mixed-gender 
studies in simulation. To our knowledge this is the first simula- 
tion study addressing the ECG signal by spectral analysis (LF/HF), 
thereby allowing a more thorough approach to stress.2° 

In the present study, no participant developed PTSD after the 
experiment. PTSD has been described in individuals exposed to 
exceptional events during emergency medicine practice with an 
immediate acute stress response.’? To our knowledge, no sim- 
ulation study has sought out an occurrence of PTSD in the SBE 
environment. We supposed that using standardized debriefing 
with ‘good-judgement method’°? might contribute to the absence 
of PTSD. Moreover, due to our Simulation Center policy regard- 
ing death no simulated patient died. Indeed, emotional reactions 
to simulated death can have a serious psychological impact, par- 
ticularly when the death results from actions or inactions by the 
learner.’! The non-occurrence of death of the simulated patient in 
our experiments may have helped to limit risk of PTSD. 


4.4, Limitations 


Despite the selection of participants and the strict protocol, this 
study had some limitations. Paramedics were older than EPs and 
PGYs. This difference should not be a bias since age does not seem 
to interfere in stress mechanisms beyond adolescence’? whereas 
level of experience should interfere in stress response.°° All par- 
ticipants were filmed and observed, which may have increased 
their stress level.’* This effect was nonetheless limited by a safe 
environment”* with a briefing explaining the presence of the 
observers. The participants were also reminded that during debrief- 
ing there was no person in judgement. On the other hand, some 
participants said they were bothered by the electrodes set up to 
record Holter parameters, which may have modified their level 
of stress. However, this factor would not have changed the rel- 
ative variation of stress parameters allowing for study of stress 
response. 


5. Conclusion 


The present high-fidelity simulation scenario tested on twelve 
multidisciplinary teams was a source of stress for all team mem- 
bers. As suspected, we found a gender difference in stress response. 
There was no correlation between STAI, Holter parameters, and 
salivary cortisol, thereby reinforcing the hypothesis of complex 
stress physiology. This study strongly favours the use of combined 
objective and subjective measures of stress to investigate stress 
response in simulation-based education. Moreover, use of such 
measures might include components of the Hypothalamic Pitu- 
itary Adrenal axis and the Autonomic Nervous System. Gender 
should be treated as an independent variable in all mixed-gender 
studies involving high-fidelity simulation. Stress is an important 
component associated with simulation-based education designed 
to approach real-life conditions in emergency situations by creating 
a realistic learning environment. Future studies should investigate 
the impact of stress response on performance. Simulation repre- 
sents a useful and effective means of promoting and advancing best 
practice in critical situations. Stress assessment could help to esti- 
mate difficult levels of tasks and training could enhance participant 
performance under stressful conditions. 
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